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^ (57) Abstract: A feedback cancellation system for a headng aid or the like adapts a first filter (1112) in the feedback path that 
^ models the quickly varying portion of the hearing aid feedback path, and adapts a second filter (1114) in the feedback path that is 
^ used either as a reference filter io£ constrained adaptation or to model more slowly varying portions of Ae feedback path. Ibe second 

filter is updated only when the hearing aid signals (1110) indicate that an accurate estimate of the feedback path can be obtained. 
Q Clianges in the second filter are then monitored (21 8) to detect changes in the hearing aidi feedback path. The first filter is ad^tively 
^ updated at least when die condition of the signal indicates that an accurate estimate of physical feedback cannot be made. It may be 

updated on a Continuous or fiiequent basis. 
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ADAPTIVF. RFFERKNrK TpR MFrHANTT^IMQ 



BACKORnTTTSfP OF TffP TTs TVFNTTnNT 
FIELD OF TH E TNVKNTTnM- 

The present invention relates to apparatus and methods for feedback cancellation 

adapted to the detection of changes in the feedback path in audio systems such as 
hearing aids. 

PRIOR ART 

Mechanical and acoustic feedback limits die maximum gain that can be achieved 
m most hearing aids. System instability caused by feedback is sometimes audible as a 
continuous high frequency tone or Whistle emanating from the hearing aid: Mechanical ■ 
vibrations from the receiver in a high power hearing aid can be reduced by combining • 

the outputs of two receivers mounted back to back so as to cancel the net mechanical 
moment; as much as 10 dB additional gain can be achieved before the onset of - 
oscillation (or NvhisUe) when this is done. But in most instruments, venting the BTE, I 

earmold or rre sheU establishes an acoustic feedback path that limits the maximum 
possible gain to less than 40 dB for a small vent and even less for large vents. The 
acoustic feedback path includes the effects of the hearing aid amplifier, receiver, and ; 
microphone as well as the vent acoustics. 

The traditional procedure for increasing the stability of a hearing aid is to reduce, 
the gain at high frequencies. Controlling feedback by modifying the system frequency 
response, however, means that the desired high frequency response of the instmment 
must be sacrificed in order to maintain stability. Phase shifters and notch filtere have 
also been tried, but have-not proven to be very effective. 

A more effective technique is feedback cancellation, in which th^ feedback 
signal is estimate and subtracted from the hucrophone signal. F^^^ .T 
typically uses an adaptive filter that models the dynamically changing feedback path 
within the hearing aid. Particularly effective feedback canceUatipn schemes are . 
disclosed in Patent Application Serial Number 08/972,265, entitled "Feedback 
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Cancellation Apparatus and Methods " incorporated herein by reference and Patent 
Application Serial Number 09/152,033 entitled "Feedback Cancellation Improvements; 
incorporated herein by reference (by the present inventors). Adaptive feedback 
cancellation systems, however, can generate a large mismatch between the feedback 
5 path and the adaptive filter modeling the feedback path when the input signal is narrow 

band or sinusoidal. Thus some adaptive feedback cancellation systems have combined 
an adaptive filter for feedback cancellation with a mechanism for reducing the hearing 
aid gain when ia periodic input signal is detected (Wyrsch, S.. and Kaelin, A., VA DSP 
implementation of a digital hearing aid with recruitinent of loudness compensation ?md 
10 acoustic echo cancellation", Proc. 1997 DEEE Workshop on Applications of Signal 

Processing io Audio and Acoustics, New Paltz, NY, Oct 1 9-22, 1 997). This approach, 
however, may reduce the hearing aid gain even if the adaptive filter is behaving 
correctiy, thus reducing the audibility of desired sounds. 

A feedback cancellation system should satisfy several performance objectives: 

1 5 The system should respond quickly to a sinusoidal input signal so that "whistiing" due 

to hearing aid instability is stopped as soon as it occurs. The system adaptation should 
be constrained so that steady state isinusoidal inputs are not canceled and audible 
processing artifacts and coloration effects are prevented from okx:urring. The system 
should be able to adapt to large changes in the feedback path that occur, for example, 

2 0 when a telephone handset is placed close to the aided ear: And the system should 

provide an indication when significant changes have occurred in the feedback path and 
are not just due to the characteristics of the input signal. 

The prefenred feedback cancellation system satisfies the above objectives. The 
system uses constrained ad^tation to limit the amount of mismatch that can occur 

2 5 between the hearing aid feedback path and the adaptive filter being used to model it. The 

constirained adaptation, however, allows a limited response to a sinusoidal signal so that 
the system can eliminate *'whistiing" when it occurs in the hearing aid. The constraints 
greatly reduce the probability that the adaptive filter will cancel a sinusoidal or narrow 
band input signal, but still allow the system to track the feedback patii changes that 

3 0 occur in daily use. The constrained adaptation uses a set of reference filter coefficients 

that describe the rnost accurate available model of the feedback path. 

Two procediires haye been developed for LMS adaptation with a constraint on 
the norm of the adaptive filter used to model the feedback path. Both approaches are 
designed to prevent the adaptive filter coefficients froni deviating too far from the 
3 5 reference coefficients. Ill the first approach, the distance of the adaptive filter 
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coefficems fron, the reference coefficients is detennined. and the nonn of adaptive 
m|er coefficient vector is clamped to prevent the distance fh>m exceeding a preset ■ 
; approach, a cost fonction is used m 

excessive deviation of the adaptive mter coefficients ftom th^ 

5 / /^'^P^^" Clamp: The feedb^^ 

adjust the FIR filter that models the feedback path (Hgu^ 

adaptation). The processing is most conveniently implemented in block time ddmain 
fonnwiththeadaptiveco^fficien^^ 

LMS adaptation adapts the mter coefficients w,(m) over the block of data to 
10.. the error signal given by 
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n=0 



where s„(m) is the microphone input signal and v„(m) is the output of the FIR filter 
modeling the feedback path for data block in. and there ar^ N samples per block The 

LMS coefficient update is given by 



1= - . (2) ■ 



w.(m + l) = w,(m) + x e„(m)g„:,(tn). 



n=0 



where g,.,(m) is the input to the adaptive filter, delayed by k samples, for block m. 

^ In general, one wants the tightest bouiid on the ad^^^^^ 

sun allows the system to adapt to expected changes in the feedback path such as those 
caused by the proximity of a telephone handset. The bound is needed to prevent 
coloration artifacts or temporary instability in the hearing aid which can often result 
from unconstrained growth of the adaptive filter coefficients in the presence of a 
sinusoidal or narrow band input signal. The measurements of the feedback path indicate 
that the path response changes by about 10 dB in magnitiide when a telephone handset 
IS placed near the aided ear. and that this relative change is independent of the type of 
earmold used. The constraint on the norm, of the adaptive filter coefficients can thus be 
expressed as , ' 
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K-l . : . 

£ Iw^(m) - w^(0)l 

lc = 0 

L iw,(0)i 

K=0 . 



1 0 



where Wj^(m) are the current filter coefficients, Wj^(O) are the filter coefficients 
determined during initialization in the hearing aid dispenser's office, the FIR filter . 
consists of K taps, and g « 2 to give the desired headroom above the reference 
condition. The clamp given by Eq (3) allows the adaptive filter coefficients to adapt 
freely when they are close to the initial values, but prevents the filter coefficients from 
growing beyond the clamp boiindary . . 

Adaptation with Cost Fimction: The cost function algorithm minimizes the error 
signal combined with a cost function based on the magnitude of the adaptive coefficient 
vector: 



N-1 



K-1 



(4) e(m) = £ [s,(m) v„(m)]^ + p £ (w,(m) - w,(0)]^ 



h=0 



k«0 



where fi is a weighting factor. The new constraint is intended to allow the feedback 
cancellation filter to freely adapt near the initial coefficients, but to penalize coefficients 
that deviate too far from the initial values. 

15 The LMS coefficient update for the cost function algorithm is given by 

;*(5)-.. 

N-1 

w,(m + 1) = wjm) - 2^p[w,(m)-w^(0)]-h2^i e„(m)g„_,(m). 

The modified LMS adaptation uses the same cross correlation operation as the 
conventional algorithni to update the coefficients, but combines the update with an 
2 0 exponential decay of the coefficients toward the initial, values. At low input signgd or 

cross correlation levels the adaptive coefficients will tend to stay in the vicinity of the 
initial values. If the magnitude of the cross correlation increases, the coefficients will 
adapt to new values that minimize the error as long as tiie magnitude of the adaptive 
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coefficients remains close to that of the initial values. However, large deviations of the 
adaptive filter coefficients from the iniUal values are prevented by the exponential decay 
which IS constantly, pushing the adaptive coefficients^^ 

nius the exponential decay gi^eatiy reduces the occurrence of processing artifacts that ' 
can result from unMnded growth in the magnimde of the ^^^^^ 

Additional references that are relevant to this appUcatibn include: 

PCt International AppUcation Publication Na 9.926.453. pubUshed in May 
; 1999; . 

PCT International Application Publication No. 9.960.822. published 
November, 1999; 

PCT international Application Publication No. 9.951.059. pubUshed October : 
1999; 

Wyrsch;SigisbertandAugustiCaelin."ADSPImplementatibnofaDigi^ v 

Hearing Aid with Recmitinent of Loudness Compensation and Acoustic 

Echo Cancellation," Workshop on Applicatibns of Signal Processing to 
Audioand Acoustics. 1997, 1-4; 

Lindemann, Eric. 'The Continuous Frequency I^amic Range Compressor." 
- IEEE Workshop on Applications of Signal Processing to Audio and 
V Accoustics,NewPaltz.NY,October 19-22, 1997; 
Czyzewski, A.. R. KrolikoWski. B, Kostek. H. Sklrzyns'ki, and A. Lorens. 

"A Method for Spectral Transposition of Speech Signal Applicable in 
Profound Hearing t^ss." IEEE Workshop on Applications of Signal 
Processing to Audio and Accoustics, New Paltz, NY October 1 9-22 

1997; , ■ '■ . 

Haykin, Simon. Adaptive Filter Theory. 3rd Ed., Prentice Hall, 1996 170- 

Kates, James M. "Feedback Cancellation in Hearing Aids: Results f^om a 

Computer Simulation," IEEE Transactions on Signal Processing 
39(3). March 1991. 553-562. . 

A need remains in the art for apparatus and methods to eliminate "whistiing" in 
unstable hearing aids while providing an accurate estimate of the feedback path. 



SUMMA R V OF THF TNyPMnr^xj 
present invention comprises a new approabh to improved feedback 
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. . cancellation in hearing aids. The approach adapts a first filter that models the quickly 
varying portion of the hearing aid feedback path, and adapts a second filter that is used 
either as a reference filter for constrained adaptation or 16 model more slowly varying 
portions of the feedback path. The first filter &at mo^^^^^ 
5 the feedback path is adaptively updated on a continuous basis. The second filter is 

updated only when the hearing aid signals indicate that an accurate estimate of the 
feedback path can be obtained. Changes in the second filter are then monitored to detect 
changes in the hearing aid feedback path. 

All audio system, such as a hearing aid, according to the present invention, 
1:0 comprises ajnicrophone or the like for providing an audio signal, feedback cancellation 

means which includes means for estimating a physical feedback signal of the audio 
systeni and means for modelling a signal processing feedback signal to conipensate for 
the estimated physical feedjback signal, an adder connected to the microphone and the 
output of the feedback cancellation for subtracting the signal processing feedback signal 
is from the audio signal to form a compensated audio signal; aiidio system processing 

means, connected to the output of the subtracting means, for processing the 
compensated audio signal, and means for estimating the condition of the audio signal 
and generating a control signal based upon the condition estiinate. The feedback 
cancellation means forms a feedback path from the output of the audio system 
2 0 processing means to the input of the subtracting means and includes a reference filter 

and a current filter^ wherein the reference filter varies only when the control signal 
. indicates that the audio signal is suitable for estimating physical feedback, and wherein 
the current filter varies at least when the control signal indicates that the signal is not 
suitable for estimadng physical feedback: 

2 5 In some embodiments, the current filter varies more frequently than the 

reference filter, usually continuously . This occurs in embodiments wherein the 
feedback signal is filtered through the current filter and the current filter is constrained 
by the reference filter. 

The current filter may only be adapted when the control signal indicates that the 

3 0 signal is not suitable for estimating physical feedback, in embodiment's wherein the 

feedback signal is filtered through the current filter and the reference filter, and the 
current filter represents a deviation applied to the reference filter. 

Frequently the means for estimating the condition of the audio signal coihprises 
means for detecting whether the signal is broadband, and the reference filter varies only 
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when the control signal indicates that the signal is broadband. For example, the audio 
system processing means computes the signal speptrum of the audio signal, the means 
for estmiati^g computes the ratio of the minimum tp the mk!^^ 

density and generates a control signal based upon the ratio,and the control si 
mdicates the audio signal is suitable when the ratio exceeds a predetermined^ 
As anothei- example, the audio system processing meahs computes the correlation : . 
matnx of the audio signal, the means for estimating computes the condition number of 

the correlation matrix and generates a control signal based upon the condition number 
and the control signal indicates the audio signal is suitable when the condition number 
falls below a predetermined threshold. / 

In the preferred embodiment, the reference filter is monitored to detect 
significantchangesinthefeedbackpathoftheaudibsystem. Also, constraining means 
prevents the current filter (or the reference filter combined with the deviation filter) from 
deviating excessively firorn the reference filter. 

. ■ ■ \ ' . ' ' ' ' • ^ • 

PRIEF DESCRIPTTON np t he PR A wmn. *? 

Figure 1 is a block diagram of the first embodiment of the present invention 
wherein the reference coefficiemA^ector is allowed to adapt under c^^^ 

Figure 2 is a flow diagram showing the process implemented by the 
embodiment of Figure 1. 

Figure 3 is a block diagram of a second embodiment of the present invention 
(simplified from the embodiment of Figure D, wherein the reference coefficient vector 
IS more simply updated by being averaged with the feedback path model coefficients. 

Figure 4 is a flow diagram showing the process implemented by the 
embodiment of Figure 3. 

Figure 5 is a block diagram of a third embodiment of the present invention 
(sirmlar to the embodiment of Figure 1, but utilizing a more parallel structure), wherein 
the reference coefficient vector is allowed to adapt under certain conditions. 

Figure 6 is a flow diagram showing the process implemented by the 
embodiment of Figure 5. 

. ■» . . . ' • 

Figure 7 is a block diagram of a fourth embodiment of the present invention 
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(simplified from the embodiment of Figure 5), wherein the reference coefficient vector 
is more simply updated by being averaged with the feedback path model coefficients. 

Figure 8 is a flow diagram showing the process implemented by the 
embodiment of Figure 7. 

■ Figure 9 is a block diagram of a fifth enibodim^ 
(similar to the embodiinent of Figure 1, but utilizing a probe signal), wherein the 
reference coefficient vector is allowed to adapt under certain conditions. 

^ ; Figure 10 is a flow diagraim showing the process implemented by the 
embbdiment of Figure 9/ 

; . Figure 11 is a simplified block diagram illustrating the basic con 
present invention. 

DESCRIPTION OF THE PREFERRED EMBODINffi>r^ 

, Figures 1 , 3, 5, 7, and 9 illustrate various embodiments of the present 
invention, while Figures 2, 4, 6, 8, and 10 illustrate the dgorithms performed by the 
embodiments. Similar reference numbers are used for similar elements between 
Figures 1, 3, 5, 7, arid 9 and between Figures 2, 4, 6, 8, and 10. 

Figure 11 is a simplified block diagram illustrating the basic concept of the 
present invention. The system includes a signal processing feedjback cancellation block 
1116 designed to cancel out the physical feedback inherent in the system. Adder 1 104 
subtracts feedback signal 11 18, representing the physical feedback of the system, from 
audio input 1 102. The result is processed by audio processing block 1 106 
(compression or the like) and the result is output signal 1 108. Audio output signal 
1 1 08 is also fed back and filtered by block 1116. 

Feedback cancellation block 1116 comprises two filters, a current filter 1112 
anjd reference filter 1 1 14. Reference filter 1 1 14 is updated only when a signal 11 10, 
indicating the condition of the audio signal, indicates that the signal condition is such 
that an accurate estimate of the feedback path can be made. Current filter 11 12 is 
updated at iezist when the signal 1 1 1 0 indicates that the audio signal is not suitable for 
an estimate of the feedback to be inade. This is the case when reference filter 1114 
represents the feedback path estimate that is made when the signal is suitable, and 
current filter 1112 represents the deviation from the more stable reference filter 1114, 
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which may be required to compensate for a sudden change in the f^^^^ 
(caused, for example, by the presence of a tone). Current filter feedback signal J 108 is 
then filtered through both cunent filter (or devi^on filter) 1112 and slower vaiying 
filter 1114 (see Figures 5 and 7). . 

Feedback cancellation, in which the feedback signal is estirria^^ 

from the nucrophone signal Js not discussed in detail herein. Feedback canceUati^^ 
OT>ically uses an adaptive filter that models the dynamically changing feedback path 

withm the hearing aid. Particularly effective feedback canceUatioh sche^^^ 
disclosed in Patent Application Serial Number 08/972.265. entiUed "Feedback 

CaiceUation Apparatus and Methods." incorporated herein by reference and Patent 
Application Serial Number 09/152.033 entitled "Feedback Cancellation Improvements " 
incorporated herein by reference. 

In oAer embodiments (see Figures 1 and 3). rrfercnce filter 1 1 14 still repres^^^^ 
the feedback path estimate that is made when the signal is suitable, but current filter 
1112 represents a frequentiy or continuously updated feedback path estimate: Feedback 
signal 1 108 is filtered only by Current filter 1112. but current filter 1 1 12is constrained 
not to deviate too drastically from reference filter 1114, 

Figure 1 is a block diagram of the first embodiment of tiie present invention 
wherein the reference coefficient vector is allowed to adapt under certain conditions' 
Figure 2 IS a flow diagram showing tfie process implement^ by the embodiment of 
Figure 1 . The improved feedback cancellation system shown in Figure 1 uses 
constrained adaptation to prevent the adaptive filter coefficients 132 fi^om deviating too 
far from the reference coefficients set at initialization. However, the reference 
coefficient vector 134 is also allowed to , adapt; it can thus move from the initial setting 
to a new set of coefficients in response to changes in the feedback patii. Coefficients 

132 used to model the feedback path adapt continuously, reacting to changes in the 
feedback path as well as to feedback "whisUing" or sinusoidal input signals. Reference * 
coefficients 134. on theotiierhand. adapt slowly or intermittently when conditions 
favorable to modeling the feedback patii are detected, and do not adapt in response to 
"whistimg" or to narrow band input signals. The reference coefficients 134 are much 
more stable dian the current feedback path model coefficients 132; the changes in 
reference coefficients 134 can tiierefore be konitored to detect significant changes in the 
feedback path such as would occur when a telephone handset is positioned close to the 
aided ear. , . 
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Figure 1 sho>ys the first embodiment of the present invention utilized in a 
conventional/hearing aid system comprising an input microphone 104, a fast Fourier, 
transform block 1 12, a hearing aid processing block 1 14, an inverse fast Fourier 
transfonn block 1 16, an amplifier 118, and a receiver 120. The actual feedback of the 
system is indicated by block 124. The sound input to the hearing aid is indicated by 
signal 102, and the sound delivered to the Wearer's ear is indibated by signal 122. 

The current (continuously updated) feedback path model consists of an adaptive 
FIR filter .132 in series with a delay 126 and a nonadaptive FIR or IIR filter 128, 
although adaptive filter 132 can be used without additional filtering stages 126, 128 or 
an adaptive: IIR filter could be used instead. Error signal 1 1 0^ e 1 (n), is the difference 
between incoming signal lOiS, s(n), and current feedback path model output signal 1 38, 
vl(n): 

The reference (intermittently updated) feedback path consists of an adaptive 
filter 134 (for example a FIR filter) in series with delay 126 and nonadaptive filter 128; 
There is a second error signal 144, e2(n), which is the difference between incoming 
signal 106 and the output 140 of reference filter 134 given by v2(n). Error signal 1 10 
is used for the IMS adaptation 130 of adaptive FIR feedback path model filter 
coefficients 132, and error signal 144 is used for the LMS adaptation 136 of the 
reference filter coefficients 134. 

The error in modeling the feedback path is given by x(n), the difference between 
the tme and the modeled FIR filter coefficients! Siqueira et al (Siqueira, M.G., Alwan. 
A., and Speece, R., "Steadystate analysis of continuous adaptation systems in hearing 
aids", Proc. 1997 IEEE Workshop on Applications of Signal Processing to Audio and 
Acoustics, New Paltz, NY, Oct 19-22, 1997) have shown that for a feedback path 
modeled by an adaptive FIR filter 

(6) E[x] = R-ip, 

where p=E[g(n)s(n)] and R=E[g(n)g'r(n)]. The error in representing model filter 
coefficients will be zero if the system input 106, s(n), and the adaptive filter input 160, 
g(n), are uncorrelated. If these two signals are correlated, however, then a bias will be 
present in the model of the feedback path. For a sinusoidal input the bias will be 
extremely , large because the expected cross correlation p will be large, and the 
cdri-elation matrix R will be singular or nearly so. Thus the inverse of the correlation 
matrix will have very large eigenvalues that will greatly amplify the non-zero cross- 
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correlation. . • 

The improved feedback canceUation is designed to update the reference 
coefficients: when the bias given by Equation (6) is expected to be smaU. and to eschew 
updating the reference coefficients when the bias is expected to be large. From Equation 
(6). the bias is expected to be large when the input signal is periodic or narrow band, 

signal conditions that wiU yield a large condition number (nitio of the largest to th^ 
smaUest eigenvalue) for the correlation matrix R. The condition number is a veiy time 
consuming quantity to calculate, but Haykih (Haykin. S.. "Adaptive Filter Theory: 3n« ; 
Edition". Prentice HallrUpper Saddle River, NJ. 1996, pp 170-171) has shown that the 
condition number for a correlation matrix is bounded by the ratio of the maximum to the 
nummum of the underlying power spectral density. Thus die ratio of the input power • 
spectral density maximum to minimum can te used to estimate the condition number . 
directly from the FFT of the input signal. 

The resulting feedback cancellation algorithm is presented in Hgure 2. . 
Referring back to Figure 1, the adaptive filter coefficients 132 for the feedback path 
model are updated for each data block. The reference filter coefficients 134 are updated 
only when the correlation matrix condition number is small, indicating favorable 

conditions for tiie adaptation. The condition number 162 is estimated from FFT 1 12 of 
the input signal 106. although other signals could be used, as well as techniques not 
based on tiie signal FFT. Jf the power spectrum minimum/maximum is large, die 
condition number is smaU and the reference coefficients are updated. If the power 
spectmm minimum/maximum is small, the condition number is large and the reference 
coefficients are not updated. Returning to Figure 2. Error signal 110 is computed in. 
step 202 and cross correlated witii model input 1 60 in step 204 (block 1 30 of Figure 1 ) 
The results of ttus cross correlation (signal 150 in Figure 1) are used to update the 
current model coefficients 132. but die amount the coefficients can change is 
constiained in step 208 as described below. 

In step 220, the signal spectrum of die incoming signal is computed (e.g. in, 
FFTblock 112ofFigure 1). Step 222 computes die min/max ratio of the spectrum to 
generate control signal 162. In step 210, error signal 144 is computed (adder 142 
subtracts signal 140 from input signal 106). Step 214 cross correlates eiror 144 with 

reference input 162 (in block 136). Step 216 updates reference coefficients 134 (via 
signals 146) if (and only iO die output from step 222 indicates that die signal is of i 
sufficiem quality, to warrant updating coefficients 1 34. Step 208 uses reference 
coefficients 134 to constrain die changes to model coefficients 132 (via signals 148). 

11, 
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Finally, step 218 tests for changes in the acoustic path (indicated by significant changes 
in reference coefficients 134). 

A mbnotonically increasing function of the power spectrum minimuni/maximum 
can be used (via control signal 162) to control the fraction of the LMS adaptive update 
that is actually used for updating reference coefficients 134 on any given data block. 
Other functions of the input signal that can be used to estimate favorable conditions for \ 
adapting the reference coefficient vector include the ratio of the maxinium of the power 
spectrum to the total ppwer in the spectrum, the maximum of the power spectrum^ the 
maxinium of the input signal time sequence,, and the average power in the input time 
sequence. Signals other than the hearing aid input 106 can also be used for estinfiating . 
favorable conditions; such signals include intermediate signals in the processing 1 14 for 
the hearing impaiitnent, the hearing output 122, and the input to the adaptive portion 
of the feedback path model 160. 

A further consideration is the level of the ambient signal at the microphone 
relative to the leyel of the signal at the microphone due to the feedback. The present 
inventor (Kates, J.M., "Feedback cancellation in hearing aids: Results from a computer 
simulation**, DEEE trans. Signal Proc, Vol. 39, pp 553-562, 1991) has shown that the 
ratio of these signal levels has a strong effect on the accuracy of the adaptive feedback 
path model. In a compression hearing aid, the lower the ambient signal level tfie higher 
the gain, resulting in a more favorable l^yel of the feedback relative to that of the 
ambient signal at the microphone and hence giving better convergence of the adaptive 
filter and a more acciurate feedback path model! , Thus the rate of ad^tation of the 
reference coefficient vector in a compression hearing aid can be increased at low input 
signal levels or equiyalently for high conipression gain values. In a hearing aid 
allowing changes iri the hearing aid gain, increasing the gain will also lead to 
improvements in the ratio of the feedback path signal relative to the ambient signal 
measured at the hearing aid microphone and hence allows more rapid adaptation of the 
reference filter. This modification of the rate of adaptation of the reference coefficient 
vector for changes in the hearing aid gain would be in addition to the algorithm shown 
in Figure 2: 

The reference coefficients 134 will be an accurate representation of the slowly 
varying feedback path characteristics. Refei-erice coefficients 134 can therefore be used 
to detect changes in the feedback path, that can in turn be used to control the hearing aid 
signal processing 1 14, Examples would be to change the hearing aid frequency 
response or compression characteristics when a telephone handset is detected, or to 
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reduce the high frequency gain of the hearing aid if a large increase in the magnitude of • 
the feedback path response were detected. Changes in the norm, in one or more 
coefficients, or in the Fourier transform of the reference coefficient vector can be used 
to identify meaningful changes in the feedback path. 

The system of Figure 1 and the associated algorithm of Figui^ 2 nearly double : 
the number of aritiimetiG operations needed for the feedback cancellation when 
compared to a system that does not adapt the reference filter coefficients. A simpler 

system (shown in Figure 3) and algorithm (shown in Figure 4) can be used if there is 
not enough processing capacity for the complete system. In the simpler system, 

reference coefficients 334 are updated by being averaged with feedback path model 
coefficients 332 rather than by using LMS adaptation. 

Let r(m) be the spectrum minimum/iimimum for data block m. track r(m^ 
a peak detector having a slow attack and a fast release time constant to give a valley 
detector, and let d(m) denote the valley detector output with 0 ^ d(m) ^ 1 . The value of 
d(m) will converge to 1 when there have been a succession of data blocks all having 
broadband power spectra; under these conditions the feedback path model will tend to : 
converge to the actual feedback path. On the other hand, d(m) will approach 0 given a ' 

narrow band or sinusoidal input signal, and will drop to a small value Whenever it 
appears that the input signal could lead to a large mismatch between the feedback path 
model and the actual feedback path. The value of d(m), or a monotonically increasing 

function of d(m), can therefore be used to control the amount of the feedback path 
model coefficients averaged with die reference coefficients to produce the new set of 
.reference, coefficients. 

The resulting system is shown in Figure 3 and the algorithm flow chart is . 
presented in Figure 4. Figure 3 is very similar to the system shown in Figure 1 . except 
that the reference coefficients 134 are not LMS adapted, which means adder 142 and 
LMS adapt block 136 can be removed. Current feedback path model 332 is updated for 
eveiry data block, and thus responds to the changes in the feedback path as well as to a 
sinusoidal input signal. For a broadband input signal 106. the reference coefficients 

334 are slowly averaged with the feedback path model coefficients (via signal 352) to 
produce the updated reference coefficients, and the averaging is slowed or stopped 
when the input signal bandwidth is reduced (controlled by signal 362). In a 

compression hearing aid, the rate of averaging can also be increased in response to 
decreases in the input signal level 106 or increases in the compression gain. In a hearing 
aid having a volume control or allowing changes in gain, die rate of averaging can be 
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increased as the gain is increased. 

Figure 4 is very similar to Figure 2, except that steps 210 (computing the 
second error signal) and 214 (cross correlating the second error signal with the 
reference input) have been removed arid block 216 (LMS adaptive reference update) has 
5 been replaced with block 416 (averaging the reference and the current model). Block 

424 has been added to low pass filter the min/max ratio of the spectrum. The output of 
step 424 controls whether the reference coefficients are averaged with the model 
coefficients. 

In the system shown in Figure 1, the first filter is the current feedback path 
10 model and represents the entire feedback path- The second filter is the reference for the 

constrained adaptation, and the second filter coefficients are updated independently 
when the data is favorable. An alternative approach is to model the feedback path with 
two adaptive filters 532, 1 34 in parallel as shown in Figure 5. The reference filter 1 34 
in this system is given by the reference coefficients (as in Figure 1), and current (or 
15 deviation) filter, 532 represents the deviation of the modeled feedback path from the 

: reference. Note that in Figure 5 and 7, the current filter (filter 1 1 12 of Figure 1 1) is 
called a deviation filter, to more clearly identify the function of the current filter in these 
embodiments. The deviation filter 532 is still adapted using constrained LMS 
adaptation; the clamp uses the distance from the zero vector instead of the distance from 
2 0 the reference coefficient vector, and the cost function approach decays the deviation 

coefficient vector towards zero instead of towards the reference coefficient vector. - 
Under ideal conditions the reference coefficients 1 34 will give the entire feedback path 
and the deviation signal 538 out of filter 532 will be zero. Deviation filter 532 is 
adapted for every block of data, aind the reference filter coefficients 534 are adaptively 

2 5 updated whenever the input data is favoi^able. In a compression hearing aid, the rate of 

adaptation of the reference filter coefficients can also be increased in response to 
decreases in the input signal leyel or increases in the coinpressioh gain. In a hearing aid 
allowing changes in the hearing aid gain, more rapid a;daptation of the reference filter 
. wouid'occur as the gain is increased. 

3 0. A somewhat different interpretation of the deviation and reference zero filters is 

that reference filter 134 represents the best estimate of the feedback path, and deviation 
filter 532 represents the deviation needed to suppress oscillation should the hearing aid 
temporarily become unstable. With this interpretation, reference filter coefficients 134 
should be updated whenever the incoming spectmm is flat, and deviation filler 
3 5 coefficients 532 should be updated whenever the incoming spectrum has a large 

14 - ■• 
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peak/vaUey ratio. The spectram minimum/maximum ratio CM therefore . 
. control the propoition of the adaptive coefficient update vectbrs used to update the 
deviation and reference coefficients for each d^ta block! An alternative would be to use ■ 
the spectrum minimum/ maximum ratio to control a switch that selects which set of 
coefRcients is updated for each data block. 

The algoritiun flow chart for die parallel filter system of Figure 5 IS presented in 
Figure 6. This flow chart is nearly identical with the flow chart of Figure 2. The only 
difference between the two algorithms is that for the paraUel system, in step 602. output; 
538 of deviation filter 532 is subtracted from 110 by adder 508, to give tiie error signal 
510. LMS update 530 cross correlates error signal 5 10 and signal 1 60 in step 604. 

Deviation filter coefficients, 532 arc then updated in step 606 (via signals 550). 
Deviation coefficient updates arc constrained in step 608. Thus, the computational 
requircments f or the paraUel system of Figure 5 will be virtually identical with tiiose for 
the system of Figure 1. ; 

1 5 In Figure 7, the alternative system of Figure 5 has been simplified in much the 

same way that die system of Figure 1 was simplified to give the system of Figure 3. A . 
portion of deviation filter coefficients 732 is added to reference filter coefficients 734 
whenever conditions are favorable. As in the case of the earlier simplified system of 
Figure 3, favorable conditions arc based on the output 562 of the valley detected 
spectrum minimum/ m^mum ratio. The value of 562, or a monotonically increasing 
function of 562. can therefore be used to control die amount of deviation coefficients. 
732 added to reference coefficients 734 to produce die new set of reference coefficients 
734. The simplified parallel systena is shown in Figure 7, and the algorithm;flow chart 
is presented in Figure 8. 
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In step 802 of Figure 8, the combined outputs of deviation filter 732 and 
reference filter 734 form signal 738. which is subti^acted f rom input 106 by adder 708 
to form errof signal 710. In step 804, LMS adapt block 730 cross correlates error 
signal 710 with model input 160: In step 806, deviation coefficients 732 are updated 
via signals 750. The amount of adaptation is c^)nstrained in step 208 filter as described 
above. Step 220 computes die signal spectrum, step 222 computes die min/max ratio, 
and step 424 low pass filters die ratio as described earlier. In step 816. if conditions 
dictate, die reference filter 734 is replaced by an averaged Version of die reference plus 
the'deviation. 

In a compression hearing aid, die rate of averaging can also be increased in 
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r^ponse to dedreases in the input signal level 1 06 or increases in the compression gain. 
In a hearing aid having a volume control or allowing changes in gain, the rate of 
averaging can be increased as the gain is increased; The computational requirements for 
this simplified system are similar to those for the system of Figure 3 since the reference 
: 5 : and deviation filter coefficients can be combined for ea^^ 
filtering operation. 

The adaptation of the reference coefficients can be improved by injecting a noise 
probe signal into the hearing aid output. Figure 9 shows the system of Figure 1 with 
the addition of a probe, signal 954. The adaptation of reference coefficients 934 uses the 

1 0 cross cbrrelatibn of the error isignal 1 44, e2(n), with the delayed, 956, and filtered, 

958, probe signal 964, g^(n). This cross correlation gives a more accurate estimate of 
the feedback path than is typically obtained by cross correlating the error signal with the 
adaptive filter input gl (n) as shown in Figure 1 . A constant amplitude probe signal can 
be used, and the adaptation of the reference filter coefficients can be performed on a 
1 . 5 continuous basis. However, a system with better accuracy will be obtained when the 

level of probe signal 954 and the rate of ads^tation of reference filter coefficients 934 
are controlled by the input signal characteristics, e.g. by signal 1 62. The preferred 
probe signal is random or pseudo-random white noise, although other signals can also 
be used. 

2 0 The probe signal amplitude and the rate of adaptation are both increased when 

the input signal has a favorable spectral shape and/or the input signal level is low. 
Under these conditions the cross correlation operation 936 will extract the maximum 
amount of information about the feedback path because the ratio of the feedback path 
signal power to the hearing aid input signal power at the microphone will be at a 

2 5 maximum. Adaptation (via signal 946) of the reference filter coefficients is slowed or 

stopped and the probe signal amplitude reduced when the input signal level is high; 
under these conditions the cross correlation is much less effectiye at producing accurate 
adaptive filter updates and it is better to hold the reiference filter coefficients at or hear 
their previous values. Other statistics from the input or other hearing aid signals as 

3 0 described for the system of Figure 1 could be used as well to control flie probe signal 

amplitude aiid the rate of adaptation. 

The adaptive Algorithm flow chart is shown in Figure 1 0. This algorithm is 
very similar to that of Figure 1^ except as follows. Gross correlation step 1014 cross 
correlates signal 964 derived from probe signal 954 with error signal 144, inLMS 
3 5 adapt block 936. In step 1016, filter 934 is updated, via signals 946. In step 1020, the 
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probe Signal level 954 is adjusted in response to the incoming sig^^ 
minimum/maximum ratio. 

What is claimed is: ... 
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1. An audio system comprising: 

nieans for providing an audio signal (1102); ' [■ . ■ 

feedback cancellation means (11 16) including means for estinriating a physical 
feedback signal of the audio system, and means for modelling a signal 
5 processing feedback signal to compensate fpr the estimated physical 

feedback signal; 

subtracting means (1 104), connected to the audio signal providing means and 
the output of the feedback cancellation means, for subtracting the signal 
processing feedback signal from the audio signal to form a compensated 
10 audio signal; 

; . aucUo system processing means (1106), connected to the 6^ 

subtrapting means, for processing the compensated audio signal; 
means for estimating the condition of the audio signal (1110) and generating a 
' control signal based upon the condition estimate; . 
15. wherein said feedback cancellation means, forms a feedback path from the output 

(1 108) of the audio system processing means to the input of the 
subtracting means and includes- 
a reference filter (1 1 14), and 
a current filter (11 12), 

2 0 wherein the reference filter varies only when the control signal indicates 

that the audio signal is suitable for estimating physical feedback, 
and wherein the current filter varies at least when the control 
signal indicates that the signal is liot suitable for estimating 
: . physical feedback, 

2. The audio system of claim 1 wherein the current filter varies more frequently 
than the reference filter.. 

3 . The audio system of claim 2 wherein the feedback signal (160) is filtered 
through the currerit filter (132); and the current filter is constrained by the reference 
filter(134). 

4 . The audio system of claim 2 wherein the current filter varies continuously: 

5 . The audio system of claim 1 wherein the feedback signal (160) is filtered 
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. through the current filter (532) and the reference filter (134); and the current filter 
represents a deviation (538) appiieid to the .reference filter. 

6. The audio systein of ciaim 1 wherein the- means for estimating the cohdition of 
the audio signal comprises means (l 12) for detecting whether the signal is broadband 

and the reference filter varies only when the control signal (162, 362. 562>i^^^^ 
that the signal is broadband. 

7 . The audio system of claim .6. Whereifl the audio system processing means 
comprises means (1 12, 220) for computing the signal spectrum of the audio signal- 
wherein the means for estimating (222) computes the ratio of the minimum to the 
maximum input power spectral density and generates a control signal (162) based upon 
the rauo; and wherein the control signal indicates the audio signal is suitable wHen the 
ratio exceeds a predetermined threshold. 

8. audio system of claim 6. wherein the audio system processing means 
comprises means for computing the correlation matiix of the audio signal; wherein the 
means for estimating (1 12. 222) computes the condition number (162) of the correlation 
matiix and generates a control signal based upon the condition number; and wherein the 
control signal indicates die audio signal is suitable when the coridition number falls 
below a predetermined threshold. 

9. The audio system of claim 1, further comprising: 

. monitoringmeans(218)formonitoringthereferencefilter(13.4)todetect . 
significant changes in die feedback path of the audio system. ' 

1 0. The audio system of claim 1 , further comprising: 

constraining means (208) for preventing the curi-em filter from deviating 

excessively.from the reference filter. ' 
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1.1. A hearing aid comprising: 

a microphone (104) for converting sound into an audio signal (106); 
feedback cancellation means (1116) including means for estimating a physical 
feedback signal of the hearing aid, and means for modelling a sighal 
5 . processing feedback signal to compensate for the estimated physical 

feedback signal- 
subtracting means (1 104)/connected to the output of the niicrophone 

output of the feedback cancellation means» for stibtracting the signal 
processing feedback signal from the audio signal to form a cpmpensateid 

1 0 audio signal; 

hearing aid processing means (1 106), connected to the output of the subtracting 

means, for processing the compensated audio signal; 
means (1 12) for estimating the condition of the aiidio signal and generating a 
control signal (1 1 10) based upon the condition estimate; and 
15 speaker means (120), connected to the output of the hearing aid processing 

. means, for converting the processed compensated audio signal into a 
sound signal; 

wherein said feedback cancellation means forms a feedback path from the output 
of the hearing aid processing means to the input of the subtracting means 
.2 0 and includes- 

a reference filter (1 1 14), and 
a current filter (1 112), 

wherein the reference filter varies only when the control signal indicates 
that the audio signal is suitable for estimating physical feedback, 

2 5 and wherein the current filter varies at least when the control 

signal indicates that the signal is not suitable for estimating 
physical feedback. 

12 . The hearing aid of claim 1 1 wherein the current filter varies more frequendy 
than the reference filter. 

13. The hearing aid of claim 12 wherein the current filter (132) represents the 
curi-ent best estimate of physical feedback; wherein the feedback signal (160) is filtered 
through the current filter; and wherein the current filter is constrained by the reference 
filter (134). . 

14. The hearing aid of claim 12 wherein the current filter varies continupusly. 
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15 . . The hearing aid of claim 11 wherein the cuircnt filter (532> r^^^ 
deviation applied to the refisrence filter; and wherein the feedback signal is filtered 
through tfie current filter and the reference filter. 

16. The hearing aid of claim 1 1 wherein the means for estimating the condition of 
the audio signal comprises meam for detecting whether the signal is broadbtod-^^^ 

reference filter varies only when the control signal indicates Aat the signal 
broadband. 

17. The hearing aid of claim 16. wherein the hearing aid processing means 
comprises means for computing the signal spectnim of the audio signal; wherein the 
means for estimating computes (220) the ratio of the maximum to minimum input 
power specttal density and generates a control signal based upon the ratio; and wherein . 

the control signal indicates the audio signal is suitable when the ratio exceeds a 
predetermined threshold. 

.18. The hearing aid of claim 16. wherein the hearing aid processing means 
comprises means for computing the correlation matrix of the audio signal; wherein the 

means for estimating computes the condition number of the correlation matrix and 
generates a control signal based upon the condition number; and wherein the control 
signal indicates the audio signal is suitable when the condition number falls below a 
predetermined threshold. 

19. The hearing aid ofclaim 11, further comprising: 

monitoring means (218) for monitoring the reference filterto detect significant 
' changes in the feedback path of the audio system. 

20, The hearing aid of clairhll .ftirther comprising: 

constiraining (208) means for preventing ihe current filter from deviating 
excessively from the reference filter. 
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